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Vortex structure of the vacuum and confinement 
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We discuss the recently published numerical computations of the vortex free energy. They dramatically demon- 
strate that the excitation probability for a sufficiently thick vortex approaches unity at large beta, i.e. that thick 
vortices 'condense' in the vacuum. This is known to imply confinement. An analytical approach for exploring 
this phenomenon is also outlined. 



1. INTRODUCTION 

A large body of recent work may be summa- 
rized as follows. Configurations carrying nonzero 
vorticity over sufficiently large scales occur in the 
vacuum with nonzero probability for all j3 < oo. 
In fact this probability approaches unity. The 
statement implies the presence of arbitrarily long 
thick vortices in the vacuum. This in turn auto- 
maticaly implies area-law for large Wilson loops, 
i.e. nonvanishing string tension for all (3 < oo. 

A vortex carrying configuration may be char- 
acterized over a sufficiently large loop by a singu- 
lar gauge transformation which is multivalued in 
SU(N), but single-valued in SU(N)/Z(N). The 
topological obstruction encountered in trying to 
extend the definition of the singular gauge trans- 
formation throughout the area enclosed by the 
loop forms then a surface of codimension 2. It is 
only the presence of the obstruction, rather than 
its exact location (which is gauge dependent) that 
is of course significant, and signals the presence 
of a 'vortex core' that cannot be deformed to the 
trivial vacuum by any continous deformation. It 
must be stressed, however, that in general such a 
configuration can have 'core' vorticity distributed 
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in a very nonlocal, patchwise manner. This in 
fact can be expected to be the case for generic 
individual configurations in the partition func- 
tion sum. Individual vortex configuration with 
spatially well-defined localized core and long dis- 
tance tail can of course be constructed and do 
occur in the path integral, but are not necessarily 
representative of most such configurations. One 
should not confuse individual vortex configura- 
tions in the partition sum with the canonical pic- 
ture of the vortex realized as a soliton of an ef- 
fective long distance action (resulting from some 
integration over configurations). 

This makes identification and computation of 
the contribution of vortex-carrying configurations 
in the path measure somewhat tricky. Thus, e.g., 
attempts to gauge fix so that vortices in any 
one configuration can be associated with well- 
localized thin 'projection vortices' encounter am- 
biguities and problems that are discussed in other 
contributions in these proceedings jjj . Here we re- 
port on recent studies H which, for the first time, 
directly measured this contribution. The idea is 
to enforce the presence of a vortex (mod N) in 
every configuration by an appropriate singular 
transformation, and measure the corresponding 
expectation. This then measures directly the con- 
tribution of such configurations, and hence the 
excitation probability (equivalently, the free en- 
ergy cost) for a vortex in the vacuum. It pro- 
vides the most direct and physically transparent 



2 



method for assessing the presence of vortices in 
the vacuum at large (3. 

2. VORTEX FREE ENERGY 

A coclosed set of plaquettes (2-cells) is a closed 
set of (d — 2)-cells on the dual lattice. Thus, in 
d = 3, it is a closed loop of dual bonds; in d = 4, 
a closed 2-dim surface of dual plaquettes. Let 
V denote a coclosed set of plaquettes that winds 
through every, say, 2-dim [12] -plane of the lat- 
tice A, i.e. a topologically nontrivial plaquette 
set wrapped around the periodic lattice (d-torus) 
in the (d — 2) directions A ^ 1, 2 perpendicular 
to fj, = 1 v = 2. We define the partition function 



Za(t) = 



(\{dU b exp( -Y.MUp) 

b ^ p£V 



pev 



(1) 



where the plaquette action A P (U P ) is replaced by 
the 'twisted' action A p {tU p ) for each plaquette of 
V. Here the 'twist' r S Z(N) is an element of the 
center. There are thus (N — 1) different nontriv- 
ial choices for r. The trivial element r — 1 is the 
ordinary partition function Z\(l) = Z^. As indi- 
cated by the notation on the l.h.s., the exact po- 
sition or shape of V is irrelevant; the only depen- 
dence is on the presence of the Z{N) flux winding 
through each [12]-plane. V can be moved around 
and distorted by a shift of integration variables, 
but not removed; it is rendered topologically sta- 
ble by winding completely around the lattice. By 
the same token introducing two twists, r on V 
and t' on V", in (|l|) is equivalent to introduc- 
ing one twist r" — tt' since V and V can be 
brought together by a shift of integration vari- 
ables. This expresses the mod N conservation 
of the Z(N) flux introduced by the twist. Thus, 
for N = 2, any odd number of such (nontrivial) 
twists is equivalent to one, and any even number 
to none. 

The twist enforces the introduction of a discon- 
tinuous (singular) SU (N) gauge transformation 
on every configuration in the path integral mea- 
sure in (|l|) with multivaluedness in Z(N); in other 
words, the introduction of a tt 1 {SU{N)/Z{N)) 



vortex. The set V represents the topological ob- 
struction to having single- valuedness everywhere. 
Thus Z\(t), eq. ([j]), is the partition sum for 
configurations with a topologically stable vortex 
trapped inside and completely winding around 
the lattice. The vortex free energy (v.f.e.) or- 
der parameter, now defined as 

exp(-F v (r)) = 
ex P ( - E { A p^ U P ) - A p {U p )) \ ), (2) 

is then the normalized expectation for the excita- 
tion of such a vortex. 

As it is evident from (|J) , measurement of this 
quantity presents a difficult sampling problem. It 
represents the expectation of an operator that 
can get significant contributions only from con- 
figurations having very small statistical weight, 
and thus becoming very difficult to measure with 
increasing lattice size. The problem can be ad- 
dressed by use of a multihistogram method. Re- 
sults of such a computation g for SU (2) are pre- 
sented in figure |l| for symmetric lattices and for 
three different values of /3. The lattice spacings 
are a = 0.165 fm, a = 0.119 fm and a = 0.085 fm 
for (3 = 2.3, (3 = 2.4, and (3 = 2.5, resp. Notice 
that, with the lattice size expressed in physical 
units, the measurements for different /3's fall on 
the same curve, as they should. This indicates 
that the universal curve has been reached, and 
will not change at larger beta. Also, the on- 
set of the sharp rise around 0.7 fm is in the re- 
gion of the finite temperature deconfining phase 
transition providing another indirect consistency 
check. The amplitude goes to unity for suffi- 
ciently large lattice size. For comparison, an up- 
per bound for Coulomb massless behavior (from 
maximizing action in spin wave approximation) 
gives - exp(-/3(7r/2) 2 ) w 0.085 at = 2.3. 

Thus the excitation free energy cost for a suffi- 
ciently thick vortex vortex goes to zero at large (3. 
The curve gives a direct estimate of 'sufficiently 
thick': slightly above I fm. This is the thickness 
at which vortices can become arbitrarily long as 
the free-energy cost goes to zero. In this sense the 
vacuuum can be viewed as having a 'condensate' 
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Figure 1. SU(2) vortex free energy (|^) vs. lattice 
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of long thick vortices. 

Over distances well below 1 fm the thick 
vortices ("confiners") are not readily detectable 
among perturbative and other shorter range fluc- 
tuations such as instantons dominating the ac- 
tion. (The instanton distribution is generally es- 
timated to be centered around 0.3 — 0.4 fm.) As a 
quantitative demonstration, let < trU p >(~> de- 
note the plaquette expectation in the presence of 
the twist, i.e. in the presence of a vortex winding 
around the lattice. For a 10 4 lattice at f3 — 2.4 
we find < trUp >(~>= 1.2557(2), while for the 
usual plaquette expectation in the absence of the 
twist < trUp >= 1.2598(2). Thus asymptotically 
in the large volume limit a vortex becomes locally 
'invisible'. 



SU (2) for simplicity. 

Consider a rectangular loop C in a [12]-plane. 
Then, for any reflection positive plaquette action, 
the Wilson W[C] loop obeys the bound f§: 



W[C] < exp(-F) 



Ac 



(3) 



where Ac is the minimal area bounded by C, A 
the area of the lattice [12]-plane, and 



exp(-F) = J2 T cxp(- J F 1 v(r)) 

T=l,-1 



l-exp(-F v (-l)) 



(4) 



defines the Z(2) Fourier transform of the v.f.e. 

Similarly, let i? m i n denote a 't Hooft loop of 
minimal length, i.e. one consisting of the 4 cubes 
(in d=4) forming the coboundary of a plaquette 
(equivalently, the four dual bonds of the bound- 
ary of a dual plaquette), and having the flux at- 
tached to the loop winding completely around the 
lattice. This is nothing but (||) with one plaque- 
tte removed from the set V on which the Z(2) 
twist resides. Physically, this may be viewed as 
an operator creating a vortex 'punctured' by a 
Z(2) monopole current loop of minimal length (a 
'tagged vortex'). Then it is not hard to prove 
that 

B min > exp(-F v (-l)) exp[-/3 < tiU p ><~>] (5) 

where, as above, < trU p >'~- ) denotes the pla- 
quette expectation in the presence of the (un- 
punctured) twist winding around the lattice. Now 
confining behavior for B mm can be shown [Q to 
provide a sufficient condition for the existence 
of a nonvanishing string tension through a lower 
bound on it. Hence, by (H), the v.f.e. does also. 



3. RELATION TO WILSON AND 'T 
HOOFT LOOPS 

The behavior of other order parameters such 
as the Wilson and 't Hooft loops can be related 
to that of the v.f.e. In particular, the existence 
of non-nonvanishing string tension is implied by 
confining behavior for the v.f.e. We consider 



4. AN ANALYTICAL APPROACH 

It would clearly be very nice to get a 
(semi) analytic handle on the v.f.e. at large /?. 
This is of course a very difficult problem. Here we 
would like to briefly outline a possible approach 
that we have been pursuing. Consider the stan- 
dard polymer expansion of the partition function. 
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It is based on the character expansion of the pla- 
quette function Boltzmann factor: 



■ trUr, 



Inserting in Z\ and expanding one obtains the 
polymer expansion: 

Za = (\{dU\{e^^ 

J p 

»(/3) iAi [En^ 

r 7 

J (/3) |A| exp(^a(C) J] 



Co ( 



Co I 



7£C 



(6) 
(7) 



The first sum (||) is over all sets T of disjoint poly- 
mers 7 of activity z 7 ; it is a finite series on a finite 
lattice A. The second sum (Q) is over all 'clus- 
ters' C, i.e. sets of overlapping polymers allowing 
multiple copies of the same polymer. There are 
simple explicit formulas for the coefficients a(C). 
This is now an infinite series even on finite A. 
As it is well-known, the free-energy cluster ex- 
pansion (Q) converges for strong coupling (small 
enough /3), and also in the |A| — > oo limit. 

The v.f.e (|^) now has a very nice expansion as 
the difference of two free energy cluster expan- 



Za(-) 
Za 



ex p ( e ac ( n 



7£C 



n 

7GC 
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Here are the activities in the presence of 

the twist. But (||) is independent of the loca- 
tion at which the twist intersects any [12]-plane. 
Hence there is a huge cancelation of terms in 
the difference: only polymers that extend com- 
pletely across [12]-planes contribute. But (||) is, 
of course, still an infinite scries even on a finite 
lattice. 

Proceeding in this manner at large (J appears 
hopeless unless one can work with a finite sum. 
Fornunately this can be accomplished by resum- 
ing all clusters with repeated multiplicities to con- 
vert the expansions (Q), (^) into a finite sum on 
a finite lattice. The result is 

ZaH 
Za 



CX P ( E E a ^c ( ^ Z^(S) - In Z(S)) ) (9) 

c ses 

where now the sum is over all skeleton clusters 
C and all sub-skeleton clusters S within each C, 
with Z(S) the partition function for S, and ex- 
plicitly computable coefficients asc- @ is an ex- 
act expression that holds for all (3. 

Having reached (||) as a suitable starting point, 
the next step is to define a decimation procedure 
for the sum over skeleton clusters on the lattice 
A in terms of skeleton clusters on a coarser lat- 
tice A'. Various decimation schemes for sets of 
clusters may be devised. If a sufficiently accu- 
rate such scheme can be implemented, it should, 
by successive iterations, reach the strong coupling 
regime without encountering a fixed point. 

5. CONCLUSIONS 

The question of confinement at large (3 can be 
reduced to that of the excitation probability for 
an (arbitrarily long) vortex. Our recent numeri- 
cal simulations have shown that this probability 
indeed tends to unity. Much work remains to be 
done. In particular, investigations on lattices of 
asymmetric length in the various directions are 
needed in order to obtain better understanding 
of the process of flux spreading leading to vortex 
thickening, as well as allow extraction of numeri- 
cal values for the string tension directly from the 
v.f.e. Related work on 't Hooft loops, free energy 
derivatives and strong coupling effective actions 
is being pusued by several groups jy . 
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